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Abstract 

The purpose of this Letter is to present a computation of the 
interaction energy of gauged 0(3) Chern-Simons vortices which are 
infinitely separated. The results will show the behaviour of the in- 
teraction energy as a function of the constant coupling the potential, 
which measures the relative strength of the matter self-coupling and 
the electromagnetic coupling. We find that vortices attract each other 
for A > 1 and repel when A < 1. When A = 1 there is a topological 
lower bound on the energy. It is possible to saturate the bound if the 
fields satisfy a set of first order partial differential equations. 
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In recent years the study of vortices has become a subject dealing with 
both particle physics and condensed matter physics. The physics of the re- 
cently discovered high-T c superconductors is an important problem. The 
studies of Chern-Simons (C-S) solitons can be related to the unusual be- 
haviour of this new type of superconductor. 

The purpose of this Letter is to present a computation of the interac- 
tion energy of gauged 0(3) Chern-Simons vortices. We discuss the solutions 
analytically and numerically, and present the numerical results. The corre- 
sponding problem for vortices of the Abelian Higgs model was first considered 
by Jacobs and Rebbi [|I|], as expected, our results are analagous to these. The 
corresponding problem for Abelian Chern-Simons Higgs vortices was consid- 
ered in [g]. 

There is an example in which this crossover, from attractive to repulsive 
behaviour occurs, namely, the case of Abelian Skyrme CP 1 vortices ||. In 
this case no Bogomol'nyi bound occurs, yet crossover takes place. 

The essential feature of C-S solitons of all 2 + 1 dimensional theories is 
that the time component of the Euler-Lagrange equation arising from the 
variation of the gauge connection is used to eliminate the time component 
of the gauge connection from the stress tensor density in the static field 
configuration. The energy, which is the integral of this static stress tensor 
density, is bounded from below by a topological charge and is finite. The 
static C-S solitons then are the solutions of the static second order equations 
of this two dimensional subsystem. The problem of gauging sigma models, 
which was considered sometime ago by Fadde'ev |3J] both in 2 + 1 and in 3 + 1 
dimensions, is a pertinent one both physically and for its own sake. Recently 
there has been renewed interest in this area, thus far in 2 + 1 dimensions. 
First the CP 1 model was gauged with a Chern-Simons £7(1) field in Refs.|J, 
and soon after that the £7(1) gauging of the CP 1 model with a Maxwell field, 
and including a Skyrme term was presented in Ref.||. 

More recently Schroers[]8| gauged the 0(3) sigma model with a Maxwell 
£7(1) field, much in the same spirit as in the earlier work of Ref.[|J. The 
0(3) sigma model was gauged with Chern-Simons dynamics in |TTJ. It is 
interesting to ask if the model supports attractive and repulsive phases. 
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1 The 0(3) Model 



We start by denning the Lagrangian on 2 + 1 dimensional Minkowski space, 

£ = ^= er'F^A, + {D,4> a f ~ (1) 

where the three component field (fi a = (4> a ,(p 3 ) with a = 1,2, is constrained 
by <p a (j) a = 1, and the covariant derivative D^(j) a is defined as in ||, and with 
the opposite sign of the coupling used in ||, as 

D^ a = d^ a + A„s a ^, D^ 3 = d^ 3 . (2) 

The Lagrangian ([!]) is U(l) gauge invariant by virtue of 2 . The potential 
function, not yet specified, is allowed only to depend only on the U(l) invari- 
ant component <f> 3 of (fi a . The energy momentum tensor of this Lagrangian 
is 

Tfj, v = 2D^ a DA a ~ g^u(D^ a D u r - W) (3) 
The Hamiltonian density is given by; 

T 00 = n = A)0 a A)0 a + A0°A0 a + w + (A <f) a ) 2 . (4) 

The essential feature of all Chern-Simons solitons of all 2 + 1 dimensional 
theories is that the time component of the Euler-Lagrange equation arising 
from the variation of the gauge connection is used to eliminate the time 
component of the gauge connection from the stress tensor density in the 
static field configuration. In the static limit this is solved for A yielding; 

K e ij Fi- 



In the static limit T 00 = H reduces to; 
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(6) 



The choice of potential function will be dictated by the requirement that the 
volume integral of 7i is bounded below by a topological charge. This fixes the 
potential uniquely. The topological charge density is not a total divergence, 



3 



but is only locally a total divergence. The usual winding number density is 
go = EijS abc di(f) a dj(j) b (j) c \ which is related to its gauged version 

The lower bound on the volume integral of the static Hamiltonian density 
(|6p can be inferred from the following inequalities; 

[ £ij D4 a - e ahc D^ h ^f > 0, {^—F l3 -V2e tJ V) 2 >Q. (7) 

V2|0 a | 

When the expansion of the static Hamiltonian density is examined, the po- 
tential is specified uniquely by the requirement that the energy density be 
bounded below by a total divergence. The lower bound can be arranged to be 
equal to the winding number density qq plus a total divergence by identifying 
the potential uniquely as; 

^=^(l-0 3 ) 3 (l + 3 ) 

The requirement that the surface integral should vanish is satisfied by choos- 
ing appropriate boundary conditions. In the case of topologically stable 
solutions, the asymptotic conditions 

lim 3 = -1, Jim 3 = ±1 

\S\— *0 >oo 

guarantee that the volume integral of g yields a nonzero integer winding 
number. This statement assumes that Ai does not grow too fast at infinity, an 
assumption which is amply justified as will be seen below when we specialise 
to the radial field configuration. The conditions (|8p with the upper sign 
pertain to the topologically stable solutions of nonzero winding, while those 
with the lower sign to the nontopological vortices. In the latter case, the 
winding number vanishes for all vorticities N, in this case the energy is given 
a lower bound by the magnetic flux. 

2 Self Dual Solutions 

The topological inequality is saturated when the inequalities (0) are sat- 
urated (for the sake of definiteness we have chosen the value of k — 1 ), 
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(8) 



yielding the Bogomol'nyi equations, 



F^ = ^-(l-0 3 ) 2 (l + 3 ) (9) 

SijD^ = ±£ a6c ^0V c , (10) 

where the lower/upper signs pertain to anti/self-duality. Our radially sym- 
metric Ansatz for the fields A, and 4> a is 

_ a(r) - N A 

— (11 j 

Q = sin /(r) n Q , 3 = cos /(r) (12) 

where Xi = — and n a = (cos iV0, sin N9) are unit vectors, with N defined to 
be an integer. 

The Bogomol'nyi equations (H) and (p70|) now reduce to the following 
pair of coupled nonlinear first order ordinary differential equations; 

^ = ±(l-cos/) 2 (l-cos/), /' = ±^£. (13) 

The topological asymptotic conditions, namely (§) with the upper sign, 
which were chosen in anticipation of our restriction to the antiselfdual case, 
now read 

lim/(r) = 7r, lim f(r) = (14) 

which for the field configuration (O) and (|T^) imply vorticity = —N. 
This is the same as in the usual (ungauged) 0(3) model where the radially 
symmetric antiselfdual vortices satisfy the asymptotic conditions ([14]) , 
while the selfdual vortices satisfy instead lim r ^ f( r ) = 0; lim^oo /(r) = ir. 

The asymptotic behaviour of the function a(r) in ( |TT|) is of no conse- 
quence to the topological stability of the soliton, unlike in the cases of the 
Higgs models 0, and of the gauged CP 1 models j|. This is because in the 
latter systems the topological charge, which is again related to the vorticity, 
is also proportional to the magnetic flux. In the case of the gauged 0(3) mod- 
els, defined in || and here, the magnetic flux of the solution is not restricted 
by the requirement of the stability of the soliton. The only constraint on the 
large r behaviour of a(r) here is the requirement that the surface integral 
bounding below the static Hamiltonian density, should vanish. This means 
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that a(r) should not grow faster than the quantity (cos / — 1) in that region. 
Since the magnetic flux is proportional to the quantity [— a(oo) + a(0)], we 
shall seek solutions for which both 0(00) and a(0) are finite, since it is rea- 
sonable that the solutions we seek correspond to finite magnetic flux field 
configurations. As explained above, we shall take a(0) = N, but will take 
a(oo) = a, where a is a nonzero constant whose sign will depend on whether 
we are considering the topological or the nontopological solutions. Corre- 
sponding to the asymptotic conditions (0) for the function f(r), we state 
the asymptotic conditions on the function a(r) as 



lim a(r) 
1 — >o 



lim air] 



a. 



(15) 



3 Non-Self dual solutions 



We will not study further the Bogomol'nyi equations, as they have been am- 
ply remarked upon in Refs. |T0| [[11]]. It should be noted that in the latter 



reference, ||11|| , the proof of existence of the solitons has been demonstrated. 
The N = 1 soliton exists, in contrast to the proof of nonexistence of the 
N = 1 soliton which has been given in || in the case of the 0(3) sigma 
model with a Maxwell term describing the curvature. 

We will examine the Euler-Lagrange equations arising from the static 
Hamiltonian. Using the Ansatz the static Hamiltonian reduces to the one 
dimensional subsystem: 



L = r 



defined by 



r sin / 



a sin / 



+ A (sin/(l-cos/)) ; 



(16) 



d 2 xH = 2tt / dr L. (17) 
When the one dimensional Lagrangian is varied the resulting equations are; 



a(r) sin 4 / + 2 cot(f)a r f r H — - — a rr 

COS f Sill f 

Aor(l — cos /) sin / (cos / — cos 2 / + sin 2 /) H 

r 

cot / esc 2 /a 2 _ _ 
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where f r = 4-. We first examine the equations fll8|) in the region r C 1 
region. The solution is; 



f(r) 
a(r) 



n + f r N + hr N+2 + 0(r N ^) 
N + A r 2N+2 + A ir 2N+4 + 0(r 




2N+6 



) where, 



(18) 
(19) 



If N 



(20) 



UN > 2 



(21) 



The constants f and A Q are fixed by the asymptotic value of the solution at 
infinity. They are found by using numerical methods, which require correct 
values of the fields at the boundary. 

Using a power series solution of the Bogomol'nyi equations, it is found 
that the ratio 



holds. This gives one free parameter in the analytic series which controls the 
behaviour of the solutions. It will also serve as a check on the numerically 
evaluated constants. 

Now considering the region r ^> 1 and anticpating decaying solutions, we 
linearise the Euler Lagrange equations about their asymptotic values. That 
is, the equations are linearised in the functions F(r) about / = and A(r) 
around the asymptotic value a(r) = a. The Euler Lagrange equation for 
f(r) is linearised about its asymptotic value of zero, and found to be; 



In order to have finite energy solutions the constant C2 is chosen as zero. The 
equation for a(r) is also linearised, and found to be; 



a 2 F - rF r - r 2 F„. = 0. 



The solutions to this are ; 



F = — + c 2 r a . 
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whose solution is a = c^r 2 or that a is a constant. The latter option is 
chosen, in order to restrict attention to finite energy solutions. This constant 
is chosen to be a. In the case of the analysis of the Bogomol'nyi equations 
there is found to be a restriction of the value of the constant a. It satisfies 
the inequality; a > 1. 

4 Numerical Solutions 

We have studied the equations numerically using a shooting method. The 
integration started in the region r<C 1 using as initial data the power series 
solutions. The constants f and A Q have been found which give the correct 
behaviour, as r ^> 1, of the functions f(r) and a(r). The profiles for the func- 
tions f(r) for vorticity N = 2,4 and those for a(r) are given in figures 1,2,4,5, 
and the respective energy density profiles in figures 3 and 6. In table 3 the 
total energies of the vortices have been calculated for the approximate solu- 
tions. If the Bogomolny'i equations are substituted into the one dimensional 
system the result is; 

2/ r asin/ 2a r 

L = 1 (1-cos/) (23) 

r r 

The total energy of a topological self dual solution can be written as; 

E sd = An J Q dr — (o(l-cos/)) (24) 
= 8n\N\ =2n(4\N\). (25) 

This last line is found using the boundary conditions for both functions a(r) 
and f(r). The energy is independent of the choice of the number a. It is 
simply enough that a r (r — > oo) = 0. The figures for the energy density, in 
both graphs and tables, are given in units of 2ir. The figures of the energy 
density for the value of constant A = 1 is the self dual limit, 4|7V| in our 
convention. It can be seen that the calculated values are close approximations 
to the analytic calculation. In figure 7, the values of total energy are plotted. 
For the n — 2 vortex, these figures are multiplied by 2. It can be seen that 
the n = 4 vortex energy is larger than the n = 2 energy, for the regime where 
Aq < 1. Also it can be seen that the n — 2 vortex energy is larger than the 
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n = 4 energy, for the regime where Ao > 1. There is a, "cross over", between 
attraction and repulsion behaviour. The relations; 



2S(n = 2, A < 1) < £{n = 4, A < 1), 2S(n = 2, A > 1) > £{n = 4, A > 1), 

can be seen from the data, and in the graph. This states that a vortex 
with degree 2 is lighter that a vortex with degree 4, so it is energetically 
favourable that vortices repel. Clearly in the second relation, attraction is 
favourable. The boundary conditions for the vortices were chosen to be the 
anti-selfdual configuration. When the vortices are self dual, that is when 
Ao = 1 the energy of 4 vortices is twice that of 2 vortices. Only at the value 
of the constant Ao = 1 are the energies per unit vortex number equal. This 
means that there is no interaction energy between self dual vortices. It can 
be shown that the diagonal components of the stress tensor vanish in the self 
dual configuration. The values of the constants A Q and f a for the self dual 



solution should be related by equation (|22[) . It is seen that the numerically 
calculated numbers are close to their theoretical ratio. 

For the nontopological solution, different boundary conditions are used 
and the result of calculating the energy is E S( i = 4n + 4a. Although the value 
of the energy in the nontopological case depends on the surface integral 
contribution and is, in fact, the magnetic flux, the energy is arbitrary. We 
have not examined the nontopological solitons, because of this. 

When the present work was completed, the work of Ref. |12] was brought 
to my attention. Our approaches to the problem of the gauged 0(3) sigma 



model are widely different, in particular, in Ref. []12| , the examination of 
varying a different a coupling constant is made. The present work is different, 
in the sense that it was carried out in the spirit of Ref. Jl| , where the Higgs 
coupling constant is varied. 
Acknowledgement : 

It is a pleasure to thank Prof. D.H. Tchrakian, who suggested the prob- 
lem, and for continued helpful discussions. 



9 



References 

[1] L. Jacobs and C. Rebbi, Phys. Rev. B19 (1979) 4486. 

[2] K. Arthur, Phys. Lett. B356 (1995) 509. 

[3] D.H. Tchrakian and K. Arthur, Phys. Lett. B352 (1995) 327. 

[4] L.D. Fadde'ev, Lett. Math. Phys. 1 (1976) 289. 

[5] M.A. Mehta, J.A. Davis and I.J.R. Aitchison, Phys. Lett. B281 (1992) 
86; 

B.M.A.G. Piette, D.H. Tchrakian and W.J. Zakrzewski, Phys. Lett 
B339 (1994) 95. 

[6] see for example A. Jaffe and C.H. Taubes, Monopoles and Vortices, 
Birkhauser, Basel, 1980. 

[7] J. Hong, Y. Kim and PY. Pac, Phys. Rev. Lett. 64 (1990) 2230; 
R. Jackiw and E. Weinberg, ibid. 64 (1990) 2234. 

[8] B.J. Schroers, Phys. Lett. B356 (1995) 291. 

[9] PK. Ghosh and S.K. Ghosh, Phys. Lett. B366 (1996) 199. 

[10] K. Kimm, K. Lee and P. Lee, preprint 3U-TP-712, hep-th 9510141. 

[11] K. Arthur, D.H. Tchrakian and Y. Yang, "TOPOLOGICAL AND 
NONTOPOLOGICAL SELFDUAL CHERN-SIMONS SOLITONS IN 
A GAUGED 0(3) SIGMA MODEL", submitted. 

[12] J. Gladikowski, preprint DTP 95-75. 



10 



Constants 


n = 


: 2 


A 


A 




0.8 


-7.90372 


-5.0000 


0.9 


-8.09182 


-5.0000 


1.0 


-8.33333 


-5.0000 


1.1 


-8.53843 


-5.0000 


1.2 


-8.65862 


-5.0000 



Constants 


n = 


-= 4 


A 


A 


F 


0.8 


-4.68536 


-5.00000 


0.9 


-4.85028 


-5.00000 


1.0 


-5.00000 


-5.00000 


1.1 


-5.13790 


-5.00000 


1.2 


-5.26629 


-5.00000 



Constants 
A 

0.8 
0.9 
1.0 
1.1 
1.2 



Energy 



E(n = 2, A) 

7.9408 

7.9694 

8.0000 

8.0411 

8.0819 



E(n = 4, A) 

15.9423 

15.9695 

16.0000 

16.0330 

16.0680 



11 



Figure 1: Profile of the function f(r) for the vortices with n = 2 with Ao = 
1.2, ...,0.8. 

Figure 2: Profile of the function a(r) for the vortices with n = 2 with lower 
values of a corresponding to lower A . 

Figure 3: Profile of the energy density for the n = 2 vortices where increasing 
peaks represent increasing energy and increasing Ao- 

Figure 4: Profile of the function /(r) for the vortices with n = 4 with Ao = 
1.2,..., 0.8. 

Figure 5: Profile of the function a(r) for the vortices with n = 4 with lower 
values of a corresponding to lower Ao- 

Figure 6: Profile of the energy density for the n = 4 vortices where increasing 
peaks represent increasing energy and increasing Ao- 

Figure 7: Graph of the energy of two superimposed vortices, 2 x £(A, n = 2) 
and £(A, n = 4), as a function of A . 
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